Massachusetts Institute of Technology 
C. S. Draper Laboratory 
Cambridge, Massachusetts 


MEMO 


TO: 

R. Larson 

FROM: 

L. Berman 

DATE: 

6 July 1970 

SUBJECT: 

Errors in RIO - Latvel and Forvel Computations 


Attached is a derivation of the exact computation of Latvel and 
Forvel, and a comparison of the approximation used in RIO. 

The conclusion drawn from this analysis is that the error in 
RIO does not lie in the resolution approximation. 




Analysis of RIO Latvel and Forvel Computation - L. Berman ' 

Current Algorithm 

The algorithm used to compute VHL (Latvel) and VIIF (Forvel) in RIO 
consists of two steps. First the moon-rate-corrected velocity, VMP, is 
resolved along horizontal axes parallel to and normal to the CSM plane: 


VHY = VMP - UHYP 
VHZ = VMP . UHZP 


( 1 ) 


where UHYP is normal to the CSM plane. 

Next, using the assumption that the platform is aligned with its y-axis 
parallel to UHYP, the velocity is transformed into body-axis related horizontal 
components by rotating thru the outer gimbal angle (o): 


VHL = VHZ sin O + VHY cos O 
VHF = VHZ cos O - VHZ sin O 


( 2 ) 


If, in fact, the platform is aligned with its y-axis along UHYP, 


UHYP = (0, 1, 0) 


(3) 


If, in addition, the alignment is at the landing site, and we are near the 
landing site, 


UHZP = (0, 0, 1) 


(4) 


Then, from (1) 


VHY = VMPY 
VHZ = VMPZ, 


(5) 


where we refer to the components in platform coordinates, and, finally, 


VHL = VMPY cos O + VMPZ sin O 
VHF = -VMPY sin O + VMPZ cos O 


( 6 ) 



! 


Correct Computation 

A correct computation depends on the definition of "forward" and "lateral". 
For example forward could be the intersection of the vehicle X-Z plane with the 
horizontal plane, or it could be the projection of the vehicle Z-axis on the 
horizontal. With both pitch and roll not zero, these two directions differ. 
Perhaps other definitions could be used. We will use the second definition in 
this analysis. 

The simplest approach is to compute a body-related set of horizontal axes 
from the body axes, since these are available, and then resolve VMP into this 
system directly: 


let UY = ZMB x UR 
and UZ = UR x UY 


(7) 


where UZ is now a horizontal projection of the body Z-axis, and UY is normal to 
it. Then 


VHL = VMP * UY 
VHF = VMP - UZ 


( 8 ) 


From pg. 5. 6-43 of R-567 (Rev 7) we can deduce 


/ sin O sin M cos I + cos O sin I 
ZNB = -sin O cos M ( 

-sin O sin M sin I + cos O cos I 


(9) 


If we again assume we are at the point of alignment: 


UR = (1, 0, 0), 


( 10 ) 


then UY = (0, -sinOsinMsinl + cosOcosI, sinOcosM) 
and 


UZ = (0, -sinO cosM, - sinO sin M sin I + cos O cosl) 


(ID 



Finally, using (8), we get 


VHL = VMPY (-sinOsinMsinl + cosOcosI) + VMPZ(sinOcosM) 

and (12) 

VHF = VMPY (-sinOcosM) + VMPZ(-sinOsinMsinl 4- cosOcosI) 

If we assume I=M = 0, 

VHL= VMPY (cosO) + VMPZ (sinO) 

(13) 

VHF=VMPY (-sinO) + VMPZ(cosO) 

which is identical with (6). 

Error in Current Algorithm 
If we assume I=M=5. 7° 

cosl = cosM = . 995 
sinl = sinM = . 100 

and from (12) 

VHL= VMPY(-. 01 sin O + . 995 cos O) + VMPZ (. 995 sinO) 

(14) 

VHF=VM PY (-. 995 sin O) + VMPZ (-.01 sin O + . 995 cos O) 

It can be seen that the error cannot get beyond about 1 or 2% of the velocity 
components. In fact, it is apparent that large inner and middle gimbal angles 
would be required to make the error significant, and even then it would be only 
a significant fraction of the actual components. 

Thus we conclude that the errors seen in RIO do not arise from the 
analytical approximations used. 



